Microbubble formation and pinch-off scaling exponent in flow-focusing devices 
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We investigate the gas jet breakup and the resulting microbubble formation in a 
microfluidic flow- focusing device using ultra high-speed imaging at 1 million frames/s. 



In recent experiments [Dollet et al, Phys. Rev. Lett. 100, 034504 (2008)] it was found 



that in the final stage of the collapse the radius of the neck scales with time with a 1/3 
power-law exponent, which suggested that gas inertia and the Bernoulli suction effect 
become important. Here, ultra high-speed imaging was used to capture the complete 
bubble contour and quantify the gas flow through the neck. It revealed that the 
resulting decrease in pressure, due to Bernoulli suction, is too low to account for an 
accelerated pinch-off. The high temporal resolution images enable us to approach 
the final moment of pinch-off to within 1 /is. We observe that the final moment of 
bubble pinch-off is characterized by a scaling exponent of 0.41 ±0.01. This exponent 
is approximately 2/5, which can be derived, based on the observation that during the 
collapse the neck becomes less slender, due to the exclusive driving through liquid 
inertia. 



1 



I. INTRODUCTION 



Liquid droplet pinch-off in ambient air or gas bubble pinch-off in ambient liquid can 
mathematically be seen as a singularity, both in space and time™ The process that leads 
to such a singularity has been widely studied in recent year^HSl anc [ i s G f major importance 
in an increasing number of medical and industrial applications. Examples of this are the 
precise formation and deposition of droplets on a substrate using inkjet technology,^ or for 
the production of medical microbubbles used in targeted drug delivery!™^ 

For the pinch-off of liquid in gas, the dynamics close to pinch-off exhibit self-similar 
behavior, which implies that the local shape of the neck is not influenced by its initial 
conditions. The radius of the neck goes to zero following a universal scaling behavior with 
ro oc r a , where r represents the time remaining until pinch-off and a the power law scaling 
exponent.^ The scaling exponent a is a signature of the physical mechanisms that drive the 
pinch-off. The formation and pinch-off of a low- viscosity liquid droplet in air is described by 
a balance between surface tension and inertia, resulting in a 2/3 scaling exponent. ^ 4 1 6 1 19 1 20 -* 

The inverted problem of the collapse of a gaseous thread in a liquid is, however, completely 
different. Initially, a simple power law was predicted based on a purely liquid inertia driven 
collapse giving rise to a 1/2 scaling exponent. 1 7 1 21 1 22 1 However, many groups report power law 
scaling exponents that are slightly larger than l/2j^l^ 4 | 2a | 24 1 I n recent work of Eggers et 
alW^ and Gekle et alf^ it was demonstrated that a coupling between the radial and axial 
length scale of the neck 10 can explain these small variations in the scaling exponent. Based 
on a slender-body calculation it is found that ot{r) = 1/2 + (— 161nr) -1 / 2 , where a slowly 
asymptotes to 1/2 when approaching pinch-off. 

In the work of Gordillo et a/P^it has been shown that gas inertia, i.e. Bernoulli suction, 
plays an important role in the bubble pinch-off. The increasing gas flow through the neck 
results in an accelerated collapse with a = l/3 . |14 | 2e | It should be noted that the smaller 
the scaling exponent a the more rapidly the radius of the neck diminishes at the instant 
of pinch-off, since the speed of collapse ro oc ar OL ~ 1 ) where the overdot denotes the time 
derivative. 

In the work of Dollet et alW^ microbubble formation in a microfluidic flow- focusing device 
was investigated. A flow-focusing device comprises two co-flowing fluids, an inner gas and 
an outer liquid phase, that are focused into a narrow channel where bubble pinch-off occurs. 



It was found that bubble formation in a square cross-sectional channel (W x H = 20 /im 
x 20 fim) showed a similar collapse behavior giving a 1/3 scaling exponent. In that paper it 
was suggested that this exponent reflects the influence of gas inertia. However, this scaling 
exponent could not be conclusively ascribed to Bernoulli suction, due to a lack of spatial 
and temporal resolution at the neck in the final stages of pinch-off. 

In this work we study the bubble formation for extremely fast bubble pinch-off in a 
microfluidic flow-focusing channel of square cross-section, using ultra high-speed imaging at 
1 Mfps. The complete spatial structure of the bubble, including its neck, was captured. This 
allowed us to not only investigate the effect of Bernoulli suction, but also the influence of 
the constituent radial and axial length scale length scales of the neck. 

Here we find that the ultimate stage of microbubble pinch-off is purely liquid inertia 
driven. In our system, the neck becomes less slender when approaching the pinch-off, giving 
rise to an exponent a = 2/5 over almost 2 decades, which is different as compared to the 
case of bubble pinch-off in the bulk as reported by Bergmann et al.^ Thoroddsen et a/.P 
and Gekle et al. among others. 

II. EXPERIMENTAL SETUP 

The experimental setup is shown in Fig. [jji. The flow-focusing device is fabricated with 
a square cross-section channel geometry, with channel width W = 60 /im and height H = 
59 /im, as depicted schematically in Fig. [TJ to ensure that the collapse occurs in the radial 
3D collapse regime onlyp3 The device was produced using rapid prototyping techniques.^ 
A homogeneous layer of negative photoresist (SU-8) is spin-coated on a silicon wafer. The 
thickness of the layer defines the channel height. A chrome mask (MESA + Institute for 
Nanotechnology, University of Twente, The Netherlands) is used in contact photolithography 
to imprint features with sizes down to 2 /im. After ultraviolet exposure a cross-linking 
reaction starts which rigidifies the photoresist that is exposed to the light. The photoresist 
that is not exposed is removed during development with isopropanol. What is left is a 
positive relief structure which can be used as a mold to imprint micron-sized channels in 
polydimethylsiloxane (PDMS) (Sylgard 184, Dow Corning). PDMS is a transparent polymer 
which is obtained by mixing two components, base and curing agent, in a 10:1 ratio in weight. 
The mixture is poured on the mold and cured in a 65°C oven for 1 hour. The PDMS slab 
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FIG. 1. (a) Schematic overview of the setup for the study of microbubble formation in microfluidic 
flow-focusing devices. A high-speed camera mounted to an inverted microscope is used to capture 
the final moment of microbubble pinch-off. Gas pressure was controlled by a pressure regulator 
connected to a sensor. The liquid flow rate was controlled by a high-precision syringe pump, (b) 
Schematic representation of a planar flow- focusing device with uniform channel height H = 59 /im 
and channel width W = 60 /im. (c) Snapshot of a high-speed recording. The outer liquid flow Qi 
forces the inner gas flow Q g to enter a narrow channel (encircled by the dashed line) in which a 
microbubble is formed. 

with imprinted microchannels is removed from the mold and then holes are punched in 
the PDMS. The PDMS slab is oxygen plasma-bonded (Harrick Plasma, Model PDC-002, 
Ithaca, NY, USA) to a glass cover plate of 1 mm thickness to close the channels. Plasma 
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bonding creates a non-reversible bond which can withstand pressures up to a few bars. | 29 | 3Q | 
The oxygen plasma turns the PDMS channel walls temporarily hydrophilic which enhances 
fluid flow and wetting of the channel walls. After closing the device, 1/16 inch outer diameter 
Telfon tubing is connected to the inlet channels, through which gas and liquid is supplied. 

Nitrogen gas is controlled by a regulator (Omega, PRG101-25) connected to a pres- 
sure sensor (Omega, DPG1000B-30G). The gas supply pressure was 12kPa. A 10% (w/w) 
solution of dishwashing liquid (Dreft, Procter & Gamble) in deionized water is flow-rate- 
controlled using a high precision syringe pump (Harvard Apparatus, PHD 2000, Holliston, 
MA, USA). The liquid, with density p = 1000 kg/m 3 , surface tension 7 = 35mN/m, and 
viscosity 77 = 1 mPa-s, wets the channel walls. The liquid surfactant solution was supplied 
at a flow rate Qi = 185 //1/min. The Reynolds number Re = pQiR/rjWH w 26, with nozzle 
radius R = WH/ (W + H) ~ 30 /im, is low enough to guarantee that the flow is laminar. 

The bubble formation process is imaged using an inverted microscope (Nikon Instruments, 
Eclipse TE2000-U, Melville, NY, USA) equipped with an extra long working distance objec- 
tive with a cover glass correction collar (Nikon Instruments, 60 x Plan Fluor ELWD N.A. 0.70 
W.D. 2.1-1.5 mm, Melville, NY, USA) and an additional 1.5x magnification lens. The sys- 
tem is operated in bright-field mode using high-intensity fiber illumination (Olympus, ILP-1, 
Zoeterwoude, The Netherlands). 

To resolve the growth of the bubble and the extremely fast bubble pinch-off at the same 
time requires a high-speed camera that is capable of recording images at a high frame rate 
and at full resolution so that the field of view is sufficient to capture the entire bubble profile 
at sufficiently high spatial resolution. These two criteria, i.e. a short interframe time (of 
the order of 1 /is) and a sufficiently large field of view means a specialized ultra high-speed 
camera is required for this task. Hence, we use the Shimadzu ultra high-speed camera 
(Shimadzu Corp., Hypervision HPV-1, Kyoto, Japan) to capture 100 consecutive images at 
a high temporal resolution of 1 Mfps (equivalent to an interframe time of 1 /is), exposure time 
of 0.5 /xs, field of view of 200 /xm x 175 /im, and with a spatial resolution of 0.68 /im/pixel. 
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FIG. 2. Time series showing the formation of a microbubble in a microfluidic flow- focusing device 
recorded at 1 Mfps. The frame number is indicated at the left of each frame. For reasons of clarity, 
the background, including the channel structure, is subtracted. A detailed image that corresponds 
to frame 86 is represented in Fig. [3^i. The camera's field of view is indicated by the dashed line in 
Fig. [IJ. The exposure time is 0.5 /is. The scale bar in the lower right corner denotes 50 fim. 
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FIG. 3. (a) Snapshot of the high-speed recording showing the formation of a microbubble corre- 
sponding to frame 86 in Fig. [2j The scale bar denotes 25 fim. (b) System of coordinates for an 
axisymmetric bubble. The shape of the gas-liquid interface r(z) is described as a function of the 
axial-coordinate z. The bubble's volume is the volume enclosed between a and b indicated on the 
z-axis. The gaseous thread forms a neck that is concave in shape with ro and r c the circumferential 
and axial radius of curvature respectively. 

III. RESULTS 

A. Extracting the collapse curves 

In Fig. [2] a time series of the formation of a microbubble is shown, where all images are 
background subtracted to improve the contrast. The first image (frame 1) shows the bubble 
almost completely blocking the narrow channel (c/. Fig. [jj). This restricts the outer liquid 
flow and the liquid starts to squeeze the gas in the radial direction forming a neck. The neck 
becomes smaller and smaller until final pinch-off, resulting in bubble detachment (frame 93) . 
The complete contour of the bubble is extracted from the recordings using image analysis 
algorithms in MATLAB (Mathworks Inc., Natick, MA, USA). In order for precise detection 
of the contour the images were resampled and bandpass filtered in the Fourier domain to 
achieve sub-pixel accuracy. The schematic of the axisymmetric shape of the bubble with the 
axis of symmetry along the z-axis is given in Fig. |3| 

In Fig. [pi a surface contour plot of the radius of the bubble r as a function of the axial 
coordinate z and the time remaining until pinch-off r is shown. The minimum radius of 
the neck r is indicated by the dashed line. In Fig. we plot r as a function of the time 
remaining until pinch-off r = t c — £, with t the time and t c the collapse time, on a linear 
scale, whereas the collapse curve is represented on a logarithmic scale in Fig. [p. The collapse 
time is defined as the moment when the neck reaches its critical radius r = and breaks. 
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FIG. 4. (a) Surface contour plot (false color) of the formation of a microbubble. The axisymmetric 
radius of the bubble is plotted as a function of the axial coordinate z and the time until pinch-off 
t — t c — t. The dashed line indicates the minimum radius of the neck ro until final collapse and 
pinch-off at the origin, (b) The time evolution of the minimum radius of the neck for three different 
experiments under the same initial conditions, (c) The logarithm of the minimum radius of the 
neck ro normalized by the nozzle radius R = 30 /xm as a function of the logarithm of the time until 
final pinch-off r, normalized by the capillary time r cap = 28 /is. The solid line represents the best 
fit to the data showing a 0.41 ±0.01 slope. The dashed lines with slope 2/5 and slope 1/3 serve 
as a guide to the eye. The vertical dotted line marks the time closest to pinch-off measured in the 
work of Dollet et a/P The error in determining the collapse time t c is visualized as the gray area. 

From the ultra high-speed imaging results it can be found that this moment occurs between 
the last frame before actual pinch-off (frame 93 in Fig. [2]) and the first frame after pinch-off 
(frame 94) . We estimate the time of collapse with sub-interframe time accuracy by assuming 
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that the collapse exhibits a power law behavior with ro oc (t c — t) a , where the exponent a 
and the collapse time t c are a priori unknown, similarly as was done in Bergmann et alW^ 
From a best fit to the data we obtain t c = 93.3 ± 1 /xs, where the maximum systematic 
error is equal to the time between two frames. Note that the error in estimating t c results 
in a deflection of the datapoints away from a straight line bounded between the curves 
log r (r ± 1 x 10 -6 s)/R indicated by the gray area in Fig. ^p. This figure also suggests that 
two different stages during bubble formation exist: in the first stage of the collapse, all data 
was found to be well approximated by a power law r /R oc (r/r cap ) a , with a = 0.29 ±0.02. 
In the final stage, when r < r cap , a scaling with exponent a = 0.41 ±0.01 is observed, 
spanning almost two decades. 



B. Liquid inertia driven pinch-off 

Approaching the singularity at pinch-off (r — » 0) , the relative importance between viscous 
forces, surface tension forces, and inertial forces are given by the Reynolds number, Weber 
number, and the capillary number. The Reynolds number, as a measure of the ratio between 
inertial forces to viscous forces, is expressed as 

Re = , (1) 

V 

with characteristic length scale r , and characteristic velocity equivalent to the radial velocity 
of the interface r . The relative importance of inertial forces with respect to surface tension 
forces is given by the Weber number 

We=^. (2) 

7 

The capillary number represents the relative importance of viscous forces to surface tension 
forces as 

Ca=^. (3) 

7 

If we now assume that r oc r a , where the experimentally determined scaling exponent is 
close to a = 2/5, it follows that Re oc t -1 / 5 , We oc t -4 / 5 , and Ca oc t -3 / 5 . This implies 
that Re, We, and Ca all diverge approaching the singularity. Accordingly, inertial forces 
must dominate both surface tension and viscous forces, hence the final stage of the collapse 
is purely liquid inertia dominated. 
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FIG. 5. The axial radius of curvature r c (squares) decreases faster compared to the circumferential 
radius of curvature ro (bullets). Hence the neck becomes less slender approaching pinch-off. Thus 
the slenderness ratio A = T c /ro, becomes smaller when approaching the pinch-off. The radii of 
curvature are normalized by the nozzle radius R — 30 /im. The time until pinch-off is normalized 
by the capillary time r cap « 28 /is. 



In Eggers et alP^ it was shown that for a liquid inertia driven collapse both the radial 
and axial length scale of the neck are important. Hence, the time evolution of the shape of 
the neck is investigated by measuring its slenderness. The slenderness ratio A is defined as 
the ratio of the axial radius of curvature to the circumferential radius of curvature of the 
neck. The larger the slenderness ratio is, the more slender the neck is. The axial radius of 
curvature is measured by locally fitting a circle with radius r c to the contour of the neck, 
whereas the circumferential radius of curvature r is equivalent to the minimum radius of the 
neck, see Fig. [3| In Fig. [5] the time evolution of the principal radii of curvature are plotted 
on a logarithmic scale for the final stage of the collapse. It is found that the axial radius of 
curvature exhibits a power law behavior r c oc r^, with (3 = 0.53. The circumferential radius 
of curvature scales as r oc r a , with a = 0.41, as was shown before (c/. Fig. ^p). 

The axial radius of curvature is found to have the more rapidly diminishing exponent 
(/3 > a), which implies that the slenderness A = r c /r oc /r a —> for r — >> 0. In other 
words, the neck profile becomes less slender approaching pinch-off, thus, both the radial and 
the axial length scales are still important. This 3D character implies that the liquid flows 
spherically inward towards the collapsing neck. Thus, it might be anticipated that, this 3D 
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collapse can be approximately described using the Rayleigh-Plesset equation for spherical 
bubble collapsd^ 

•• 3 -2 1 ( 2 l\ 

r oro + -r = - [p , (4) 

with capillary pressure p. It should be noted that a necessary condition for this is that the 
neck should be much smaller than the channel dimensions (tq <C W, H). By substituting 
r oc r a in above equation and keeping the right-hand side constant, which assumes an 
inertia-dominated flow, it is found that it is necessary that a = 2/5, which agrees surprisingly 
well with our experimental findings of 0.41 ± 0.01. 



C. "Filling effect" 

How to account for the scaling r oc t- - 29±0 - 02 for r > r cap , i.e. at early times? At this 
initial stage of the collapse a thin layer of liquid with a thickness of several micrometers 
separates the bubble from the hydrophilic channel wall. 32 The liquid flow in such a confined 
channel can be described using Darcy's law for pressure driven flow through porous media. 
The volumetric flow rate of liquid that permeates into the neck region is 

kAdp 

Qin = (5) 

T] OZ 

with k the permeability, A = WH — nr^ the cross-sectional area of the thin liquid layer 
surrounding the bubble, and dp/dz the pressure gradient. 

The pressure distribution in the liquid is inhomogeneous, thus the bubble's surface does 
not have a constant curvature even though the gas pressure is practically uniform. The 
pressure gradient that drives the liquid flow can be derived from the capillary pressure 

P - 7 (- ~ - V (6) 



In the initial stage of the collapse, i.e. at the onset of neck formation, r c > r . As a 
gross simplification, we approximate the neck as a radially collapsing cylinder, of length r c 
much larger than its radius r . The capillary pressure is then p « 7/r , therefore dp/dz « 
— rQ 2 dr G /dz. 

The volumetric gas flow rate that is pushed out of the neck region is 

Qout = ~V g w -r r c r , (7) 
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with V g ~ r c rQ the volume occupied by the gas. 

The gas in the neck is replaced by the liquid. This is referred to as the "filling effect" , 

thus, from a balance between Eq. ^ and Eq. ([7]), we now get 

1 <9r 1 r . 

?o oz r r c 

hence, assuming that r c varies little in this initial stage of the collapse, t^Tq is roughly 
constant. It follows that the radius of the neck must scale as r oc r a ) with a = 1/3. This is 
in good agreement with the experimentally measured scaling exponent a ~ 0.29 ±0.02 for 

7" ^ 7~cap* 



IV. DISCUSSION 

In Gekle et a/P^ a supersonic air flow through the neck is visualized using smoke particles 
and it is reported that Bernoulli suction accelerates the collapse. An accelerated collapse 
due to Bernoulli suction is also reported by Gordillo et a/., 8 giving rise to a 1/3 scaling 
exponent. It is extremely difficult to measure the gas velocity in a microfluidic flow- focusing 
device in a direct way. However, the camera's wide field of view (200 fim x 175 /im) enabled 
us to capture the contour of the expanding bubble in great detail and allows for an estimate 
of the gas velocity. 

The volume of the bubble V^, as the gas volume downstream of the neck that is enclosed 
by the bubble contours, is calculated as follows: 

V h = f b dzirr 2 (z), (9) 

J a 

with the profile of the bubble r(z), with a the axial coordinate of the location of the neck and 
b the tip of the bubble (cf. Fig. [3]). We plot the bubble volume 14 as a function of time until 
pinch-off in Fig. [6j The bubble's contour is indicated for four characteristic moments during 
the bubble formation process in the panels (i-iv) of Fig. [6| In the initial stage the gaseous 
thread in front of the flow-focusing channel is forced to enter the channel and completely 
fills it (i-ii) . The volume of the bubble increases until it reaches a maximum volume of 38 pi 
at r = 46 /is (ii) . The restricted liquid flow starts to "squeeze" the gaseous thread and a 
clearly visible neck begins to develop. Then a remarkable event takes place — the gas flow 
reverses and the neck starts to collapse (ii-iii) until bubble pinch-off occurs. The volume of 
the bubble beyond pinch-off (r < 0) is 31 pi, which is equivalent to a bubble radius of 19 fim. 
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FIG. 6. Volume of the bubble Vb as a function of the time until pinch-off r. The volume of 
the bubble was calculated by integration over its contour along the z-axis between the neck and 
the tip of the bubble (indicated a and b in Fig. [3]). The insets (i-iv) depict the contours that 
enclose the bubble's volume corresponding to the marked data points in the graph. A second order 
polynomial fit is used to calculate dV^/dr (dashed line). The bubble reaches its maximum volume 
at r rev (dashed-dotted line) and the gas flow direction reverses, consequently, the bubble shrinks 
during the final moments before pinch-off. Different symbols represent different experiments, giving 
an indication of the reproducibility. 



We estimate the volumetric gas flow rate Q g through the neck as the time derivative 
of the volume of the bubble, Q g = 14, where it is assumed that no gas diffusion into 
the surrounding liquid take place. The bubble volume is approximated by a second order 
polynomial function, as indicated by the dashed line in Fig. |6j which is used to obtain the 
time derivative of the volume. 

The gas velocity through the neck u g is calculated as the volume flow rate Q g divided by 
the cross-sectional area of the neck (nr^). In Fig. [7] we plot the gas velocity as a function of 
the time until pinch-off. Note that the gas velocity is low during almost the entire collapse 
process, i.e. \u g \ < 0.5m/s (see the inset in Fig. [7]), however, in the final stage of the collapse, 
a strong acceleration of the gas flow is observed, reaching a velocity up to u max = — 23m/s. 
The pressure drop due to the fast gas flow through the neck is given by Bernoulli's equation 
as Ap —p g (u: 



j g v^rev 



^max) /2 ~ 0.3 kPa, with ^ rev = Om/s the gas velocity at the moment 
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FIG. 7. Gas velocity u g in the neck as a function of the time until collapse r. The initial positive 
gas velocity reverses its flow direction and accelerates when approaching the pinch-off. At the 
moment of pinch-off maximum velocity of —23 m/s is reached. In the inset an enlarged section 
of the graph for the data points encircled by the dashed line is represented demonstrating the gas 
flow reversal. Again, different symbols/colors represent different individual experiments. 

of flow reversal and p g = 1.2kg/m 3 the gas density. 

We now compare this pressure drop with the capillary pressure in the neck. Just before 
pinch-off, the capillary pressure, as a consequence of surface tension forces acting on the 
curved interface, should be at a much higher pressure than the surrounding liquid. For an 
axisymmetric surface profile, with r = r(z), the capillary pressure, as a function of the axial 
coordinate z, is given by the Laplace equation 



p(z) = 7 



1 



"y/l + r 



/2 



(10) 



(l + r' 2 ) 3/2 . 

where prime denotes the derivative with respect to z 33 Note that, at the location where 
the neck is thinnest, the first term equals the circumferential curvature (r^ 1 ), whereas the 
second term represents the axial curvature (r" 1 ). In Fig. [8j both the Bernoulli pressure 
drop and the capillary pressure in the neck are plotted as a function of the time remaining 
until pinch-off. The capillary pressure (represented by the bullets) is obtained by inserting 



the complex shape of the neck into Eq. (10). The dashed line and the dashed-dotted line 



indicate the pressure contribution from the circumferential curvature (p oc 7T a ) and the 
axial curvature (p oc — 7T _/3 ) respectively, while the solid line represents the sum of both 
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FIG. 8. Evolution of the capillary pressure and the Bernoulli pressure drop during the final moments 
of bubble pinch-off. The dashed line indicates the capillary pressure contribution due to the 
circumferential curvature (p = 7/ro, with tq oc t 41 ); the dashed-dotted line shows the pressure 
contribution from the axial curvature (p = — 7/r c , with r c oc r 53 ); the solid line represents the 
sum of both contributions. The bullets represents the capillary pressure obtained from the local 



shape of the neck (r(z)) extracted from the wide-field-of-view images and using Eq. (10). The 
increasing gas velocity through the neck causes a local pressure reduction in the neck, referred to 
as Bernoulli suction (p = —pUg/2), as indicated by the crosses. It can be seen that the capillary 
pressure clearly dominates during the entire bubble formation process. 

contributions of the capillary pressure. In the figure it is demonstrated that the pressure 
drop due to Bernoulli suction is marginal in comparison to the increasing capillary pressure 
approaching pinch-off. Hence, it can be concluded that Bernoulli suction, i.e. gas inertia, 
is irrelevant during the entire bubble formation process. It is also shown that the concave 
axial curvature counteracts the circumferential curvature leading to a significant decrease in 
capillary pressure. This confirms that the axial length scale of the neck is important and 
gives the collapse a three-dimensional character. 



V. CONCLUSION 



In conclusion, we visualized the complete microbubble formation and extremely fast 
bubble pinch-off in a microscopically narrow flow-focusing channel of square cross-section 
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(W x H = 60 /xm x 60 /xm) , using ultra high-speed imaging. The camera's wide field of 
view enabled visualization of all the features of bubble formation, including the two prin- 
cipal radii of curvature of the bubble's neck. Recording was performed at 1 Mfps, thereby, 
approaching the moment of pinch-off to within 1 /xs. It was found that the neck's axial 
length scale decreases faster than the radial one, ensuring that the neck becomes less and 
less slender, collapsing spherically towards a point sink. We describe this collapse using 
the Rayleigh-Plesset equation for spherical bubble collapse p2 and recover a 2/5 power law 
exponent which is consistent with our experimental findings. The gas velocity through the 
neck is calculated from the growth-rate of the bubble. Just before pinch-off the gas velocity 
accelerates up to —23 m/s reducing the bubble's volume, however this velocity is too low 
for Bernoulli suction to be the dominant effect. Thus, the final moment of microbubble 
pinch-off in a flow-focusing system is purely liquid inertia driven. 
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